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One intriguing fingerprint of surface states in topological insulators is the Aharonov-Bohm effect in mag¬ 
netoconductivity of nanowires. We show that surface states in nanowires of Dirac materials (bismuth, bismuth 
antimony, and lead tin chalcogenides) being in non-topological phase, exhibit the same effect as amendment 
to magnetoconductivity of the bulk states. We consider a simple model of a cylindrical nanowire, which is 
described by the 3D Dirac equation with a general T-invariant boundary condition. The boundary condition 
is determined by a single phenomenological parameter whose sign defines topological-like and non-topological 
surface states. The non-topological surface states emerge outside the gap. In longitudinal magnetic field B 
they lead to Aharonov-Bohm amendment for the density of states and correspondingly for conductivity of the 
nanowire. The phase of these magnetooscillations increases with B from tt to 27r. 


1. Introduction. 

In Dirac materials electron and hole energy spectra 
are described by a standard or modified Dirac equation 
(for review, see, for example [1]). Recently introduced 
concept of topological characterization for bulk energy 
structure of crystals brings additional attention to such 
Dirac materials as bismuth, bismuth antimony alloys 
and lead tin chalcogenides mia. For definite percentage 
of antimony the alloy Bii_a;Sba; becomes the topologi¬ 
cal insulator (TI), although pure bismuth is not the TI. 
The latter materials for certain content of tin HEiiniiz]) 
are known as topological crystalline insulators (TCIs). 
Surface states (SSs) in TIs and TCIs are protected by 
the time-inversion or mirror symmetry. The topological 
SSs result in the Aharonov-Bohm effect in magnetocon¬ 
ductivity of nanowire. In this paper we reveal that the 
Aharonov-Bohm effect can also exist in non-topological 
phase of the Dirac materials. 

Lead tin chalcogenides have the rocksalt crystal 
structure with four inequivalent L-valleys in the recip¬ 
rocal space. Bismuth antimony is of quite similar struc¬ 
ture with an additional shear along a spatial diagonal of 
the cube, which reduces the number of L-valleys to the 
three ones. In the L-point of the Brillouin zone the both 
materials have two closely lying spin-degenerate energy 
bands in the vicinity of Fermi energy Therefore 

minimal k-p-theory of the materials should include the 
two interacting bands. It leads to an anisotropic 3D 
Dirac equation that describes dynamics of fermions in 
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the k-p-approximation. In this paper we will use a sim¬ 
plified isotropic version of the equation 

{rriTz 0 tTo + CTx 0 (crfc)} vj/ = E'il, (I) 

where k = {kx,ky,kz) is the 3D momentum, 2m is the 
bulk energy gap between conduction and valence bands 
in the L-valley, cr = {ax,cFy, <Jz) and t = {Tx,Ty,Tz) are 
vectors of the Pauli matrices acting in spin and band 
subspaces correspondingly, ctq, tq are 2x2 identity ma¬ 
trices, c is an effective speed of light, Ik = ('I'c'I'^) is 
a four-component envelope function, Ike, are spinors 
describing conduction and valence bands respectively. 

Study of SSs for the 3D Dirac equation © has a 
rather long history. First, Jackiw and Rebbi [10] re¬ 
vealed that the Dirac equation has doubly degenerate 
bound states with zero energy if the mass term m has 
opposite sign from two sides of an interface region. Later 
it was shown m that the interface states with 2D Weyl 
dispersion should arise in lead tin chalcogenides junc¬ 
tions where variation of tin content results in inversion 
of the mass term. One may use this treatment to phe¬ 
nomenologically consider TI [H] or TCI [13] surface. 

Another approach to SSs of the Dirac equation © 
is to derive an appropriate boundary condition (BC) 
for the four-component wave function T* at a surface. 
This was done in Refs. m before the famous Ref. m, 
however the former paper has not been of much inter¬ 
est. The BC resulted from the Hermiticity of the Dirac 
Hamiltonian and time-reversal invariance of the prob¬ 
lem in restricted area of space. Such a BC reads as 
follows 

- iao(Tn^c)reS = 0 ) ( 2 ) 


1 
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where n = n(S) is an inner normal to the surface S, ao 
is a real phenomenological parameter which character¬ 
izes atomic surface structure and bulk band properties. 
Afterwords m it was shown that to simulate the BC 
© via some interface potentials one should consider not 
only variation of the mass term but an additional elec¬ 
tric potential which is responsible for electronic affinity 
of stacked materials. 

For special surface orientation we may additionally 
require the BC (I2|) to be invariant under mirror reflec¬ 
tion if there is no surface reconstruction m- For ex¬ 
ample, for L-valley, [111] edge, of the rocksalt structure, 
in coordinate system 2 ;||[ 111 ], 2 /||[ 110 ], and a:||[112], xz- 
plane is a plane of mirror symmetry. If surface normal 
is in the ccz-plane, i.e. n = (nx,0,nz), then invariance 
of the BC (2]) under the mirror symmetry leads to the 
condition 


MyT-^ 


= D(My)rD-\My)D{My)^ 


= -F^' = 0, 

s 

(3) 

where the mirror symmetry operator My is rep¬ 
resented by the matrix D{My) = —ixz 0 Uy, 
vh = D{My)^{x, —y, z), the BC operator is 


F = 


-iao{nz;ax+nzCrz) 1 
0 0 


(4) 


Comparing o, ® with the BC m one can see that 
the mirror symmetry does not impose any restrictions 
on the value of the parameter oq. 

In a halfspace z > 0 Eq. o with the BCs @ re¬ 
sult in appearing of the 2D massless SSs with a conical 
spectrum m (see insets in Figd]) 


E = sv\k\i 


En 


vm Eq 

— - S -|fc|| 

m 


> 0 , 


(5) 


where u = 2aoc/(l + Oq) is the SS speed, 
fc|| = {kx,ky,0), Eq = m{l - al)/{l -bag) is the energy 
of the Dirac point, here s = ±1 are eigenvalues of a 
chirality operator Tz ® (cr,[n,p]). Depending on the 
oo sign one can distinguish two qualitatively different 
cases: (i) oq > 0 the SSs are in the bulk gap (red lines 
in inset in FiglT}]), and (ii) oq < 0 the SSs are outside 
the gap (red lines in inset Fig®). Hence by means of 
the BC ([2]) one may phenomenologically describe the 
two types of SSs in lead tin chalcogenides existing in 
direct and inverse band order [niiig. 

Aim of our paper to study non-topological and 
topological-like SSs, for the sake of comparison, and to 
demonstrate their contribution in Aharonov-Bohm os¬ 
cillations of the nanowire conductivity. We will suppose 
that oo is constant on the wire surface. Analogous ap¬ 
proach was used to explain emergence and conducting 


nature of edge states circulating around a nanohole in 
graphene [18]. Besides we calculate contribution of the 
SSs to the magnetoconductivity. It is of direct connec¬ 
tion with recent experiments in magnetoresistance of 
bismuth [TQ] [20] and tin chalcogenide [21] nanowires. 

The paper is organized as follows. First, we calculate 
spectra of electrons in the nanowire without magnetic 
field. Second, we add a longitudinal magnetic field and 
also calculate contribution of the SSs to the magneto¬ 
conductivity. 

2. Dirac Fermions Spectra in Nanowire. 

Consider a cylindrical nanowire with radius R. We 
choose z-axis along the nanowire axis. The spinors 
'kc'kj; obey the 3D Dirac equation © and satisfy the 
BC © with the constant oq at a cylindrical surface 
of the nanowire. From now on we set h = c = 1 ev¬ 
erywhere, except where it is needed. Cylindrical sym¬ 
metry of the 3D Dirac equation o with BCs m im¬ 
plies conservation of longitudinal momentum kz and to¬ 
tal angular momentum Jz = <Jq ® jz, with eigenval¬ 
ues j = ±1/2,±3/2,..., where jz = ao{-ide) +(Jzl2. 
Hence, one can find as follows: 


/ V'ci(r-)e*(J \ . 


( 6 ) 


By use of the Dirac equation o, it is convenient 
to express via Tc. The radial wave functions 
'4’ci{r),'ijjc2{r) obey the Bessel equation 


- yf; ± V'cl,c2 = _ /j2^ 

(7) 

and BCs 

/ i (^dr - - ao{E + m)^ \ f ipci{r) \ 

y -kz i (dr + - ao{E + m)^ j ^ i’c2{r) J 

(8) ~ 

Therefore ipci,c 2 {'<’) = Ci^ 2 Jj^i/ 2 ikr), where 

k = ^E'^ — m? — kl, Jj^i/ 2 {kr) are the Bessel 

functions of the first kind, Ci , C 2 are arbitrary con¬ 
stants. The BCs m impose a relationship for the 
constants Ci , C 2 and give the dispersion equation 


k 


Jj-i/2{kR) 
. J3 + l/2{kR) 


Jj+i/2{kR) 

Jj-i/2{kR) 


= [ao - — ] E + m(ao + — ) . 

V aoj V “0/ 

(9) 

Imagine values of k correspond to the SSs, real values 
of k describe quantum-confined states. Because of a 
symmetry property of the dispersion equation m ao 
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l/ao,E —)> —E, it is enough to consider the case |ao| < 
1. Besides, the dispersion equation does not depend on 
j sign as follows from the property of the Bessel func¬ 
tion of integer index J_j^i/ 2 {x) = 

It reflects the axial symmetry of the nanowire. Conse¬ 
quently all ID subbands have double degeneracy due to 
the j sign. Result of qualitative graphic solution of the 
dispersion equation Q is shown on Fig.(IT]). The en¬ 
ergy of quantum-confined states at |A:^| —>■ oo aspires to 
E = ±hc^nn?/T9c^ + kl+ where 7j±i/2,n 

is the n-th zero of the Bessel function Jjj-i/2- At oq > 0 
the SSs acquire a mass equals hv\j\lR in the thick wire 
limit (|ao|mi?/hc ^ 1). The energy spectrum of the SSs 
has asymptotes 

E = svh^Jkl + p+Eo, 

K=^- s^x/kl + 4 > 0 

c^h m y z — 


when their decay length k~^ is much smaller than the 
radius R {kR » j^). The closer energy of the SSs to the 
region of quantum-confined states (this region is filled 
with grey color on the FiglT]) the more delocalized are 
the wave functions of the SSs. 

3. Dirac Electrons in nanowire in longitudinal 
magnetic field: Spectra and Conductivity. 

In magnetic field along the wire axis B = {0,0, B) 
we make the Peierls substitution p ^ p + eA in the 
3D Dirac equation (HD, -e is the electron charge. For 
the vector-potential we choose the cylindrical gauge 
A = {—Byf2, Bxj2, 0). Further, to be specific we con¬ 
sider the case B > 0. The radial components '0ci,'0c2 
of the spinor 'I'c obey the following equation 


d 

df'2 j-Qj. 


r2 ^ A2 



- 771,2 - kl)ycl,c2, 


( 11 ) 

where = 1/eB is the magnetic length squared. Nor¬ 
malizable solutions of the radial equations (ED are ex¬ 
pressed in terms of the Rummer’s function 
m- While B > 0 we are interested in spectra of SSs 
with j < —1/2. The dispersion equation for those j is: 


Numerical solution of the equation (fT^ yields the en¬ 
ergy spectrum of the Dirac fermions in the nanowire 
in the longitudinal magnetic field, see Fig. HI In two 
limiting cases the dispersion law of the SSs is: 

Ek^js = + Eo, 

(14) 

where $ = TreBR?/he is the number of the magnetic 
flux quanta through the wire cross section, with 7 _b = 0 
in a quasiclassical limit \kR/ j\ ^ max |/|, 4) (weak mag¬ 
netic field and strong localization for radial motion of 
the SSs), and 7 b = 1/2 in the limit of strong magnetic 
fields ($ \j — 1/2|, $ X^k"^/2, the last condition 
results in |ao -I- l/ao| 1 in the Ea. (fT4l) l. For oq > 0 
both values of s = ±1 are valid in Eqs. ED and surface 
subbands are formed in the bulk gap. In this case sit¬ 
uation resembles that of topological insulator nanowire 
EUHSlEn]- In the case oq < 0 only one value of s is 
allowed and surface subbands coexist with Landau and 
skipping orbit subbands (see FigHD). 

While the magnetic field is relatively weak, the gap¬ 
less SS subbands periodically emerge at half-integer val¬ 
ues of Jo = —4* with unit period. With further in¬ 
crease of the magnetic field when the magnetic length 
compared with the SS decay length k~^, the period 
is of weakly dependence on B. Far from Dirac point 
{E — Eo)R/hv = pR 1 

Density of SSs with specta ED is described by the 
formula: 

p{E)=po (^l + 2^S{E)^ X 

x0 [(E - F;o)sgn(ao(ao - 1)) -f hvke] , (15) 

where 0[...] is the Heaviside step function which deter¬ 
mines energy region where SSs exist ((for 0 < oq < 1 
the region is E < Eq + hvke, for oq > 1 the region 
is E > Eq — hvke, for —1 < oq < 0 the region is E > 
Eq — hvke, and for oq < 1 the region is E < Eo + hvke)), 
po = {E — Eo)/2Tr{hv)'^ is density of states (DOS) of free 
massless 2D DFs, ke = 2|oo|77i/?ic|l —agl, S is oscillating 
amendment to free DOS due to quantization of perimet¬ 
ric motion of SSs: 


2{j -1/2)-aoR{E + m)M 
R^k^ aoR{E + m) 


L2(j-1/2) 


-b 


M 


-b klR^ = 0. 


( 12 ) 


where 

~ M(l-A2fcV2,-j+3/2,R2/2A2) 

M(-A2fc2/2,-j-bl/2,R2/2A2) ' 


, , ^cos(27rpiin 

S{E) = Y^ ^ 

n—1 


j) cos(27rn ($ 
y/2npRn 




( 16 ) 

When energy in (1161) coincides with a bottom of every 
surface subband the quantities pR± (4) —1/2-b7 b) tend 
to integers N±. In this case the series dUD can be esti¬ 
mated as (see my- 
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Fig.l. Subband spectra (E = ER/hc) of Dirac fermions in a nanowire for two types of boundary parameter oq: a) at 
ao = —0.1 and mR/hc = 2 b) at ao = 0.4 and mR/hc = 3; bold red curves emphasize spectra of the non-topological 
surface subbands in Fig[IJi and topological-like ones in Fig[T}D. Spectra of SSs for half-space problem are shown on insets 
M- Grey colour fills regions of quantum-confined states. All subbands have double degeneracy due to the sign of j. 



Fig.2. Energy spectrum of Dirac fermions in a nanowire with a longitudinal magnetic field for two types of boundary 
parameter oq: a,c) Dimensionless spectra E = ER/hc of ID subband bottoms as a function of total angular momentum 
j for ao = 0.8 (a) and ao = —0.2 (c). b,d) Spectra of ID subbands for ao = 0.8, j = [—29.5, —21.5] (b) and ao = —0.2 
j = [—51.5, —7.5] (d). On Fig. [2}i red curves correspond to the non-topological surface subbands, blue curves - Landau 
and skipping orbit subbands. Numerical solution was carried out at hc/X = 0.896m, 4>/<l?o = 35. Dashed-dotted straight 
line on hgures b) and d) indicates the Fermi-level Ji. 


S{fi) 


0 (pR ± ($ - i - 7 b) - N±) 
2^2ttpR (pi? ± ($ - i - 7 b) - N±) ’ 


(17) 


It is seen that in this case S has square root singularities 
as it should be for DOS of one-dimensional subbands. 

4. Conductivity of Surface States. Let us calcu¬ 
late the electric conductivity cr of the .^-length nanowire 
in static uniform electric field of strength F. At first 
we consider the case oq > 0 when Fermi energy p is in 
the bulk gap (see FigI5)D). In this case only the SSs con¬ 
tribute to the conductivity at low temperatures. We ne¬ 
glect by intervalley scattering and consider weak scalar 
disorder potential as follows V (r) = with 


u{r — Ti) = US{r — ri)6{9 — 9i)6{z — z^/vi. We will em¬ 
ploy wave functions of SSs in the quasiclassical limit: 



s) 





/ 

g*0-i/2)e 

\ 




J® Jkl + 


^K{r— R)-\-ikzZ 







Jkl + (2^e^U-i/2)e 
kz ^ 


\/2'KK,r 


1 


J 



(18) 


Here we suppose that p lies well above Eq (but in the 
bulk gap) so that a lot of surface subbands with s = 1 to 
be hlled {QR ^ 1 with hvQ = p — Eq). In case of weak 
impurities while the characteristic time t of change of 
the electron’s distribution function much greater than 
the period of the evolution operator oscillation (^ h/p) 
we can restrict ourselves by r-approximation in the 
Boltzmann kinetic equation. We also suppose that one¬ 
dimensional localization length is much greater than the 
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nanowire length. It is easy to meet this condition while 
the number of filled surface subbands is great {QR » 1). 
Therefore one may assume the specific ID localization 
effects to be small. In the given approximation the cur¬ 
rent in z direction is expressed by the formula [22) 



1/ 



(19) 


some cumbersome integration we come to the final re¬ 
sult in leading approximation in 

cr = cr2D27ri? [ci — C2S{fj,)] (24) 


where ci ~ 0.95, C2 ~ 10.6, and 


CT2D 


e^v^hR 

irnimpU'^K 


(25) 


here u = {k^^j}, (vz)^ is the diagonal matrix element 
of the velocity operator of the Dirac Hamiltonian, is 
the number of non-equivalent L-valleys in the Brillouin 
zone, Tn is relaxation time in the state \ v) that expressed 
as follows 

, 20 , 

In the above formula transition probability W^v' aver¬ 
aged over configuration of scatterers is expressed by the 
following way 


= 

— 2ridri 

h Jo 




2 

X 


8{E,-E,,), 

where we introduce the Fourier coefficient of u{r 


( 21 ) 

- d) 


«a(?’) = ^tS- 

Ti), 


21^1 i-ll2 
^ 6{r — 


2'Kir. 


A = {m,q} . 

( 22 ) 


Using the Poisson summation formula one can obtain 
the following equation for relaxation time (001) in the 
leading order in {pR)~^J^: 




1 + 


2k'iR'^ 


1 + 


4U2(l+lfl53 


^ 2-V2+-^ 


S{E) 


(23) 


where riim = N/2'KRf is 2D concentration of impurities, 
j = j -I- $ — 7 b. Only impurities, that are located in 
layer of k~^ thickness around the surface of nanowire, 
effectively scatter the Dirac fermions in SSs. Additional 
phase shifting im in the last term in (1231) is manifesta¬ 
tion of the Berry phase of the SSs. Under our consid¬ 
eration {pR ^ 1) the sum in the above equation is a 
small oscillating amendment to non-quantized SS den¬ 
sity of states po- It is the amendment that leads to 
the Aharonov-Bohm oscillations of observable quanti¬ 
ties. To obtain conductivity of the nanowire we substi¬ 
tute relaxation time from Eq. (E!?)) into Eq. (ITII)) . After 


is conventional conductivity of 2D massless DFs with 
short-range impurities with corresponding replacement 
of potential strength U'^k/R —?> Uq [53]. Series in 
Eq. (l24l) converges for every value of magnetic flux $ but 
the fluxes when Fermi level coincides with bottoms of 
the surface subbands. As it was mentioned in the previ¬ 
ous section this occurs at integers QR± ($ — 1/2 — 7 b). 

In this case the r-approximation is not valid even for 
weak impurities. Eg. (1241) is valid only when the Fermi 
level is quite far from surface subband bottoms where 
1 . 

For the case oq < 0 similar calculation can be per¬ 
formed taking only transitions between the SSs. This 
contribution to the conductivity is of the form: 

Sa = tT2D27ri? [ci — C25'(p)] 0 [{p — Eq) sgn(l — Og) -I- Hvke] ■ 

(26) 

Heaviside step function in above equation means that 
this contribution is not zero when Fermi level lies in 
conduction or valence band, depending on the value of 
oo, and intersects surface subbands. 

Now, we discuss situation when the Fermi level lies 
in the vicinity of the Dirac point Eg at oq > 0. Two gap¬ 
less subbands arise at half integers $. Backscattering 
between the gapless states moving in opposite directions 
is forbidden as is zero in Eq. (EB). Therefore, in this 
regime conductance of the nanowire will $o-periodically 
oscillate with maximum value /h at half when mass¬ 
less subband arise |5^. (in weak magnetic field at half 
integers $, but in strong magnetic field limit at integers 
$ (see Ea. (|TT)) U 

Conclusions. We reveal that the T-invariant 
boundary condition for 3D Dirac equation in nanowire 
can lead to emergence of topological-like or non- 
topological types of SSs depending on sign of the BC 
parameter ag. This parameter is determined by micro¬ 
scopic surface structure and bulk band properties. At 
oq > 0 topological-like SSs forming in the bulk gap as 
in topological insulator nanowires [531 Hi IMIISO] may 
lead to the Aharonov-Bohm oscillations of magnetocon¬ 
ductivity. Although at oq < 0 non-topological surface 
subbands emerge out of the gap, they also may lead to 
the Aharonov-Bohm oscillations of the nanowire con¬ 
ductivity. The phase of magnetooscillations is tt for 
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both types of surface subbands (i.e. both signs of oq) 
in quasiclassical limit, but it acquires additional phase 
shift TT in strong magnetic field limit. Thus, observa¬ 
tion of the Aharonov-Bohm conductivity magnetooscil¬ 
lations (or surface contribution in zero B) in the Dirac 
nanowires is not evidence for topological nature of the 
SSs. 
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